Nonlinear stability of homogeneous models in 
Newtonian cosmology 

Gerhard Rein 
Mathematisches Institut 
der Universitat Munchen 

Theresienstr. 39 
80333 Munchen, Germany 

Abstract 

We consider the Vlasov-Poisson system in a cosmological setting as 
studied in jn| and prove nonlinear stability of homogeneous solutions 
against small, spatially periodic perturbations in the L°°-norm of the 
spatial mass density. This result is connected with the question of how 
large scale structures such as galaxies have evolved out of the homogeneous 
state of the early universe. 

1 Introduction 

In textbooks on astrophysics and cosmology the formation of galaxies is some- 
times explained by demonstrating that for certain model equations spatially 
homogeneous solutions, which are to represent the early state of the universe, 
are unstable against small perturbations so that such perturbations grow and 
eventually lead to the formation of large scale structures, cf. O, O, pl| . 

We intend to study the stability problem of such homogeneous states in the 
context of the Vlasov-Poisson system: 

d t f + v-d x f-d x U-d v f = Q, (1.1) 

AU = iwp, (1.2) 

p(t,x) = / f(t,x,v)dv. (1.3) 



Here t>0 denotes time, xGlR 3 position, and velR 3 velocity, / = f(t,x,v) de- 
notes the particle distribution function on phase space, p — p(t,x) denotes the 
spatial mass density generated by /, and U = U(t,x) denotes the gravitational 
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potential. Therefore, in this model the particles move according to Newton's 
equations of motion 

x = v, v= —d x U(t,x) 

under the influence of the gravitational field which they create collectively ac- 
cording to Newton's law of gravity, and collisions as well as relativistic effects 
are neglected. 

The corresponding initial-value problem where the phase space density / is 
prescribed at t — Q is well understood for the case of an isolated system, i. e., if 



the Vlasov-Poisson system (1.1), (|l.2|), (1-3) is supplemented with the boundary 



condition lim^^oo U(t,x) = and a corresponding boundary condition for /, cf. 
[|[ [lC], |l3|, |l9). Here we are interested in solutions which model the evolution 
of a Newtonian universe, and for this situation these boundary conditions are 
not appropriate. In the context of general relativity an isolated system can 
be thought of as a localized deviation from flat space — an asymptotically flat 
spacetime — whereas a cosmological solution is a deviation from a homogeneous 
state, usually one with non-zero, spatially constant mass density. Both situa- 
tions have been investigated with the Vlasov equation as matter model under 
the assumption of special symmetries, cf. |l6|, [lTj and the references therein. 
However, the existence theory for this so-called Vlasov-Einstein system is not 
yet sufficiently developed to form a basis for the stability analysis intended here. 

Spatially homogeneous, cosmological solutions for the Vlasov-Poisson system 
can be constructed as follows: For a nonnegative, compactly supported function 
HeC^M) we set 

f (t,x,v):=H{a 2 (t)\v-^x\ 2 ^ (1.4) 
where a is a positive, scalar function to be determined later. Then 

p (t,x) = J fo(t,x,v)dv = a~ 3 (t) J H(v 2 )dv, 
and, after normalizing 

J H{v 2 )dv = l, 
we obtain the homogeneous mass density 

p (t) = a- 3 {t), i>0. (1.5) 



A solution of the Poisson equation (L2) is then given by 

U (t,x):=—a~ 3 (t)x 2 , t>0, xeJR 3 , (1.6) 

o 

and it remains to determine the function a in such a way that fo satisfies the 



Vlasov equation (1.1) with force term 

4tt 



d x U (t,x) = ^-a- 3 {t)x. 
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A short computation shows that this is the case if a is a solution of the differential 
equation 

Air 

&+—a- 2 = 0, (1.7) 

which is the equation of radial motion in the gravitational field of a point mass. 
It is the purpose of the present paper to investigate the nonlinear stability of 
such homogeneous models (fo,po,Uo) against small perturbations. In order to 
have a meaningful concept of stability we restrict ourselves to the case where 
this homogeneous solution exists — and expands — for all time t > Q, which is the 
case if a(0) > and 

„ 1 . , ,o 4-7T , , . 

E a :=-a(t) 2 - — a- 1 (t)>0, 



i. e., the energy, a conserved quantity along solutions of (1.7), is positive. If E a = 
the homogeneous solution also expands for all future time but not sufficiently 
fast for our stability proof to go through. We normalize the initial condition by 
requiring that o(0) = 1; the energy then is positive iff 



d(0) > y/8n/3. 

Note that the condition a(0) = 1 means that the universe has already expanded 
out of the initial big bang singularity a = for some time. By the assumption 
on the energy there exist constants < c± < C2 such that 

ci<a(t)<c 2 , t>0. 

In order to study the stability of such a homogeneous model we investigate the 
time evolution of small dev iat ions f rom it, that is, we consider solutions of the 
Vlasov-Poisson system (1.1), fll.j), (1.3) of the form 

f = f +g, p=p +cr, U = U + W. 

To investigate the perturbations g, a, W, which we want to assume spatially 
periodic, it is useful to perform the following transformation of variables: 

x = a~ (t)x, v = v — a~ 1 (t)a(t)x,g(t,x,v)=g(t,x,v). (1-8) 

If we compute the system satisfied by g, er, W in these transformed variables and 
afterwards drop the tildas we obtain the following version of the Vlasov-Poisson 
system which governs the time evolution of deviations from the homogeneous 
state: 

d t g+-v-d x g--(d x W + av) ■d v g = 2aH'(a 2 v 2 )vd x W, (1.9) 
a a V / 

AW = 4wa 2 a, (1.10) 
a(t,x)— / g(t,x,v)dv. (1-H) 
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It should be kept in mind that now x,v <E IR 3 are related to the original variables 
by the transformation ( |1.8| ). More details on this transformation are given in 
[ p"8| , its necessity comes from the fact that in the original variables the equations 
for the deviation g would become explicitly x-dependent which would exclude 
the possibility of studying spatially periodic deviations, a class which seems 
physically reasonable and is convenient for our mathematical analysis. 

The main result of the present paper is that the homogeneous solutions 
described above are nonlinearly stable against spatially periodic perturbations 
in the following sense; for explanation of our notation and the set of admissable 
data T> we refer to the next section: 

Theorem 1.1 For every t\ > there exists eo > such that every solution g of 



(|L9|), ( 1.10 ), ( 1.11 ) with g<ET> and ||g||oo<£o satisfies the estimate 

|k(i)||oo<a- 3 (i)ei, <>0. 



// we define the density contrast S(t,x) of the corresponding solution of (1.1). 
(0), (|J by 

p(t,x) = p (t)(l + 6(t,x)) 

then 

||«(t)||oo<ei, t>0. 

Note that the result is really one on Lyapunov stability and should not be 
mistaken as an asymptotic stability result: the factor a~ 3 (t) in the estimate for 
er simply reflects the overall expansion of the solution; the model contains no 
dispersive effects which would make the solution / converge to /o- 

We briefly explain how the paper proceeds and how it relates to other work 
in this field: An essential prerequisite for a stability analysis is a global-in-time 
exis te nce th e orem for the corresponding initial value problem for the system 



(L£), ( 1.10 ), ( 1.11 ). This was established in 18 . For easier reference we collect 
some results and notation from that paper in the next section. Stability and in- 
stability results for the Vlasov-Poisson system have been established for various, 
geometrically different settings in the plasma physics case, cf. gl|, lH. 
In the more difficult stellar dynamics case linearized stability of isolated sys- 
tems is investigated in || , and the question of nonlinear stability was attacked 
in ^o| with dubious result, cf. also f2^| . Our paper differs from these in that 
we obtain a stability estimate with respect to the L°°-norm of p and obtain 
a rigorous, nonlinear result in the stellar dynamics case. On the other hand, 
our techniques are more akin to small data results known for the Vlasov-Poisson 
and related systems such as the Vlasov-Maxwell and Vlasov-Einstein systems in 
the case of an isolated configuration, cf. |l], |TJ, 0. The essential new difficulty 
arises from the fact that as opposed to these results we cannot use the disper- 
sive effect of the almost free streaming characteristics for small fields in empty 
space, since we are in a spatially periodic situation with a non-zero background 
field. Instead, we make use of the expansion of the background solution which 
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is — for the case considered here — fast enough to prevent any conglomeration 
and growth of small perturbations. A related result, but for a different matter 
model, namely dust, is given in ||. However, note that there the authors have 
to assume a positive cosmological constant to have a sufficiently rapid expan- 
sion of the background solution to prevent growth of the perturbations. As a 
first step towards our result we investigate in the third section how solutions, 
in particular their force field, depend on the initial data for g, that is, we prove 
continuous dependence on the initial data in an appropriate sense. Next we 
investigate how the characteristics behave under a small perturbation of the 
background field and establish certain resulting fall-off estimates for the spatial 
mass density of the perturbation. In the fifth section we combine these two sets 



of estimates in a bootstrap argument to prove Theorem 1.1: By making the 



initial data small and using the continuous dependence of the field on the ini- 
tial data we make sure that the field is a small perturbation of the background 
field in the sense of Section 4 so that by the results of that section we get new 
fall-off estimates on the perturbed field which turn out to be asymptotically 
stronger than the smallness assumption we start with. The time interval on 
which these fall-off estimates hold can therefore be extended to all of [0,oo[, and 
the proof is complete. Finally, we will briefly discuss how our result relates to 
the usual textbook explanation of the formation of galaxies. If this explanation 
is mathematically false, then why do galaxies exist? 



2 Notation and preliminary results 

Let Q := [0, l] 3 and S :=Q x 1R 3 . We need the following spaces of periodic func- 
tions: 

V{Q) := {h-.Wt? -^JR\h(x + a) = h(x), zGlR 3 , aeZ 3 J, 
V{S) := (h:M 6 ^TEL\h(x + a,v)=h(x,v), x,v&M 3 , aeZ 3 }, 

C"(Q) :=C"(iR 3 )nP(g), 

:= C n (lR 6 )nP{S), 
Q C (S) := {h£CZ(S)\3u>0:h(x,v) = 0, |v|>u}. 

For pS[l,oo] we denote by || • || p the usual L p -norm, where the integral (or 
supremum) extends over Q or S as the case may be, and 

L?(Q):=D>(Q)nV(Q), I%(S) :=Zf{S)nV(S), 

where we identify functions in V(Q) and V(S) with their restrictions to Q or S 
respectively. For easier reference we collect some results from |ll| : 
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Theorem 2.1 For every initial datum g£C^ v (S) with J s g — and g(x,v) + 
H(v 2 )>0 for x,vgTR 3 there exi sts a u ni que solution g S C 1 ([0,oo[x]R 6 ) (~l 
C([0,oo[;Cl(S)) of the system Q, Cl0|) , (plj ) u«tfi g(0)=g. 

This is Thm. 6.1 of [^8[ ; note that in the present case the homogeneous back- 
ground solution exists for all t>0. The restriction J s g=0 means that the 
perturbation does not change the total mass of the background solution and 
is necessary in order to allow for spatially periodic solutions of the Poisson 
equation ( l.lOj ), g{x,v)+H(v 2 ) > means that /(0) is a nonnegative function. 

Next we have the following representation of g in terms of the characteristics 
of (|L~9|), cf. Lemma 4.3 and Lemma 4.4]: 

Lemma 2.2 Let (X,V)(s,t,x,v) denote the solution of the characteristic system 
1 



l-(d x W(s lX ) + d( S )v), (2.1) 
lis V / 



a(s) ' a(s) 
with (X ,V)(t,t,x,v) = (x,v) , x^efR 3 , 0<s<t. Then 

g(t,x,v)=g((X,V){0,t,x,v))+H(V 2 (0,t,x,v))-H(a 2 (t)v 2 ), t>0, x^efR 3 . 

Finally, we need the following estimates on the solutions of the Poisson equation 
in the spatially periodic set-up: 

Lemma 2.3 For a gC^(Q) with Jga(x)dx = there exists in C 2 (Q) a unique 
solution of 

AW = Ana 

with jQW(x)dx — 0. It satisfies the following estimates: 

IIWU^C^HlpH^ + Hr), (2.2) 

II^IU < c((H-||a|| 00 )(l + lii*||Mloo/2) + |k||i), ( 2 -3) 
||9>|| 0O <C'(|| C r|| 1 /16 || f 7|^ 4 ||9 !E <7||^ 16 + ||a|| 1 ), (2.4) 



where 



In* s : 



s , s < 1 
1 + lns , s>l 



Except for the estimate (2.4) this is contained in |lj, Lemma 4.1 and Lemma 
4.2], the estimate (2.4) follows from jj], Lemma 1] and the fact that the singu- 
larity in the Green's function for the present, spatially periodic situation is the 
same as for the problem on the whole space, cf. |l8|, Lemma 4.1]. 
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In what follows we consider only small perturbations g which will be taken 
from the set 

V:=^geCl c (S) | J g = 0, g(x,v) + H(v 2 )>0 for all x,v<ElR 3 , 
||ff||i,oo <do, and g(x,v) = for |v| >u 
Here d >0 and w >0 are some fixed constants, and ||<?||i,oo denotes the L° 



norm of g and its first order derivatives. Note that in Theorem 1.1 we need to 
make ||g||oo small but we also need a uniform bound on || <7 1| i,oo ■ I n order to 
keep our notation simple we use the following 

Convention on constants: Constants denoted by C may only depend on 
the initial data set V, on c\ and C2 determined by < c\ < a(t) <C2< oo, t > 0, 
and, starting from Section 4, also on the parameter S introduced there, and 
they may change from line to line. The same convention applies to functions 
CeC([0,oo[). 



3 Dependence of the field on the initial data 

The purpose of this section is to show that we can make the perturbed field and 
its spatial derivative small on any given, finite time interval by making 1 1 <7 1 1 oo 
small. As a first step in this direction we establish this result for the perturbed 
density and field: 

Lemma 3.1 The re e xis ts a no ndecr easing function C £C([0,oo[) such that for 
every solution of (|l.9| ), ( 1.10|) , ( 1.11 ) with initial condition g£T> the estimate 

Mt)\\i + |k(*)||oo + ||0*W(t)||oo < C(t) ||.g||oo 

holds for t>0. 

Proof: There exists a nondecreasing function Pe C([0,oo[) such that 

g{t,x,v) = Q, \v\>P{t), rrelR 3 , t>0 (3.1) 

for every solution with initial condition g£T>. This follows from | fL8| , Thm. 6.1]; 
that this function can be chosen uniformly on the set V follows from the fact 
that all parameters which enter into its construction are uniformly bounded on 
T>. Lemma 2.1 and the mean value theorem imply that 

||s(t)||ao<]|s||ao+C sup \V(0,t,x,v)-a(t)v\. 

x,v£TR 3 



Since 

d_ 

ds 



(a(s)V(s,t,x,v)^ = -d x W(s,X(s,t,x,v)) 
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it follows that 



||<?(*)l|oo<||<?||oo + C / \\d x W (Snoods. 



With (3.1) we conclude that 



and 



<j(t)\\ 00 <C(t)[\\g\\ 00 + / \\d x W(s) \\cods 



RW(s)||oo<is 



\Ht)h<C(t) (\\g\\ 



(3.2) 



(3.3) 



Inserting these estimates into (2.2) yields the Gronwall type inequality 



RW(*)lloo <C(t) I || fl °||oo + J \\d x W(s) 
from which the estimate 



follows. If we insert this into ( |3.2| ) and ( |3.3| ) the proof is complete. □ 

Next we want to prove an analogous estimate for the x-derivatives of the per- 
turbed field. To do so we need certain estimates on the cc-derivatives of the 
characteristics which will follow from the following slightly rewritten form of 
the characteristic system (2.1): If 

y(s) :=a(s)X(s,t,x,v) 

for t,x,v fixed then y satisfies the differential equation 



y- 



47T 1 1 

Y^) y ~^) dxW ^ y/a{s ^ 



(3.4) 



where we have used (1.7) to express d. 



Lemma 3.2 There exists a constant C such that every solution of (2J) satisfies 
the estimates 

\d x X(s,t,x,v)\ < C^4exp (cj o- x (t) WdlW^Udr 



\d x V(s,t,x,v)\ < 



a(s) 



||SgW r (r)|| 00 |a s jr(T,t,i,i;)|dT 



for x,ve1R and 0<s<t. 
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Proof: Let £(s) = a(s)d x X(s,t,x,v) = d x y(s). Then by (3.4) we have 
and 

f(t) = a(i) id, f(i) = a(t)id. 

Thus 
and 

<a(t) + <*,(*-*) + / / (Ca- 3 (r) + a- 2 (r) 119^)11^)1^)1 drda 
<Ca(t)+J^ j ^(Co^W + o-'WHflgWWHooJI^TjIdT 

< Ca(t)+J^ r(Ca- 3 (r) + a-V)l|9xW^)lloo) |C(r)|dr. 
Now a Gronwall argument yields 

< Ca(t)expfc7 / ra' 3 (r) dA exp ( [ tcl' 2 {r)\\dlW {r^dr 



which is the first estimate of the lemma. Next observe that 

V(s,t,x,v) = ^-v+-^- f d x W(r,X(r,t,x,v))dr; (3.5) 



a(s) a(s) 
this is the solution of the initial value problem 

y^^^lv ^d x W(s,X(s,t,x,v)), V(t) = v. 

a(s) a(s) 



1 '« 



Therefore, 

d x V(s,t,x,v) = ^— I d$.W(r,X(T,t,x,v)) d x X(r,t,x,v)dT 
from which the second estimate of the lemma follows. □ 
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Lemma 3.3 There exists a nondecreasing function CgC([0,oo[) such that for 
every solution of (1.9), ( l.lOD , (1-H) with initial condition g£V the estimate 

\\d 2 x W(t)\\oo<C(t)\\g\\^ 16 

holds for t>0. 



Proof: Lemma 2.2 implies that for g£~D, 



\\d x g(t)\\oo < do sup l\d x X(0,t,x,v)\ + \d x V(Q,t,x,v) 
+ C sup \d x V(0,t,x,v)\. 

x.veTR 3 



From Lemma 3.2 it follows that 



\d x X(0,t,x,v)\<Ca(t)exp [C / \\d x W {T)\\ x dT 



and 



\d x V(0,t,x,v)\ <cj^ W^W^Woo^expfcJ a- x (r) \\d 2 x W (r^d^j ds 
< Ca(t) f liagW^JHooOcpfc / ||flgW(T)||oo^ ds 



Therefore, 



< Ca(i)exp C / Wd^W^W^dT . 



||^)]]oo<C(*)exp C/ WdlWiT^dT 



and by ( p.l|) , 



||0*<r(*)lloo<C(t)exp [C / II^M^MIIoodr 



If we insert this estimate into (2.3) and observe that ||cr(£)||i and ||cr(t) || oo can be 
bounded by some C(t), cf. Lemma 3.1, we obtain the Gronwall type inequality 



C(i)exp [C / WdlW^Udr 



\\%W(t)\\oo <C(i)+C(i)ln* 



<c(t)+c?(t) / Wdlw^Udr. 

Jo 

A Gronwall argument implies that there exists a function CgC([0,oo[) such 
that for solutions with g £ V the estimate 

|[W*)||oo <<?(*), *>0, 
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and thus also 



\\d x a(t) IU<C(t), t>0, 



holds. If we now insert this estimate and the estimates from Lemma 3.1 into 
(2.4) the assertion of the lemma follows. □ 



4 The free streaming condition and its conse- 
quences 

Although this terminology is not quite correct we refer in the following to the 
characteristics of the homogeneous background solution as the free streaming 
characteristics. We need to make precise what it means that a solution behaves 
almost like the homogeneous background solution. To this end we fix a constant 
<5g]0,13/16[. For T>0 and a parameter 7g]0,l] we say that a solution g of 



([L9), (1.10), (1.11) satisfies the free streaming condition on the time interval 



[0,T] with parameter 7>0 if 

||$«W(t)||co<7a~ 5 (*) 
(FS 7 ) { , 0<i<T. 

\\%w(t)\\oo<-ra- s (t) 

We want to study the behaviour of the characteristics of (11.91) under the assump- 



tion (FS7). As a first step we estimate the derivatives considered in Lemma ^2 
under the assumption (FS7): 

Lemma 4.1 Let the condition (FS7) hold for some 7G]0,1] on some time in- 
terval [0,T]. Then 

\d x X(s,t,x,v)\<C^\, \d x V{s,t,x,v)\<G-i^r, 0<s<t<T, x, v eM 3 . 
a[s) a(s) 

Proof: By (FS7) and since a(r) > 1 + C]T and 7< 1, 



a- 1 {T)\\d 2 x W{T)\\ 00 dT< / a' 1 - 5 (r)dT< 00. 







Inserting this estimate into the first estimate in Lemma 3.2 proves the first 



assertion, and inserting the latter into the second estimate in Lemma 3.2 yields 



\d x V(s,t,x,v)\<^-C f\ a - s (r)^\dr<C^\ 1 [°° a^^dr, 
a(s) J s a(T) a{s) J 

and the proof is complete. □ 

Remark: Without the normalizing condition a(0) = 1 the constant C in 



Lemma 4.1 would depend on a(0) in such a way that C — >oo as a(0)— >0. We 



will come back to this observation at the end of the last section. 
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Lemma 4.2 Let the condition (FS7) hold for 7G]0,1] on some time interval 
[0,T]. Then the estimate 



a(t) 



d v V(0,t,x,v)— id 



<C 7 



holds for x,weIR 3 and £6 [0,T]. J/7S]0,1] is small enough then the mapping 



is one-to-one, and 

Proof: Let 

Then 



]R 3 3w^^(0,£,x,i;)eIR 3 



|det9„V r (0,£,a;,'i;)| > -a 3 (£). 



£(s) :=a(s)9„X(s,£,a;,i?). 
4tt 1 1 



^--oT-tWs,^))^ 



3 a 3 (s) s a 2 (s) 



where i(s) =X(s,t,x,v), cf. (3.4). In the free streaming case d 2 W = the solu- 
tion takes the form 



£ F {s) = -a(s)a(t) 

this is the solution of the initial value problem 

4tt 1 



dr 



-id; 



Define 

Then 

and 

Therefore, 



3 a 3 (s) 

C(*)=C(*)=o, 



IC00I < Ca-^^K^I + Ca^a- 1 -*^) 

< Ca- 2 - 5 (s)|C(s)|+Ca(£)a- 2 - 5/2 (s) / a" 1 - 5 / 2 ^)^ 

< Ca- 2 - s (s)\((s)\+Ca(t)a- 2 - s/2 (s) 
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and 



\((s)\<Ca(t) a- 2 - s / 2 (a)da + C a' 2 ' 5 ^) \({t)\cIt, 



\({s)\ < Ca(t) / a- 2 - 5 / 2 (T)dTdcr + C / a- 2 - 5 (T)\((T)\dTdcr 

J S J (7 J S <J <7 

= Ca(t) [ [ daa~ 2 - s/2 (T)dT + C [ [ daa' 2 ^ '(t)|C(t)|<£t 



s J s 

t 



<Ca{t)J raT 1 (r)dr + C J ra-^ fl (r)|C(r)|dr 

<Ca(t) + C j a _1 - a (T)|C(r)|dT, 
from which we conclude that 

\((s)\<Ca(t)exp(c J (T^Cr) drj <Ca(t), 



1. e., 



i „ , N i aft) , , f dr 

\d v X(s,t,x,v)\<C-^r+a(t) 
a{s) 



From (p.5|) we obtain 



„ . . a(t) 1 

d v V{S,t,X,V) = — r T ld+— — 

a(s) a{s) 



a 2 (r)- 

dlW(T,x{T))d v X{T,t,x,v)dT, 



and thus 



\d v V(0,t,x,v)-a(t)id\ < 7 / a~ d (T) C-^4-+a(t) 



,a(t) 



da 



a(r) J T a 2 (a) 



<Cja{t)[ I a~ 1 ' s {r)dT+ I flT*(r) 



a 2 (a) 



dr 



< C 7 a(i) (l + J a- 1 - 6 ' 2 ^) J a-^^i^dadi^j 

< C 7 a(t), 



1. e., 



9„V(0,t,a;,u)— id 



an estimate which holds for all te [0,T], XjVEJR 3 , and 7 g]0,1]. By making 7 
small and using the mean value theorem we conclude that the mapping v * 
V(0,t,x,v) is one-to-one and 



det——d v V(0,t,x,v) 
a(t) 



> 



1 
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for all x, v G R 3 and te[0,T}. □ 

Using the estimates of the previous lemma we can now show that the spatial 
mass density of the perturbation decays in time: 



Lemma 4.3 Let 7 be as in Lemma and let g be a .solution with g(zT>, which 
satisfies (FS7) on some time interval [0,T]. Then the estimates 



and 



lk(t)lloo + lk(«)||i<Ca- 3 (*)(|| S || 00 + 7 ) 
||^a(t)|| oo <C , a- 2 (t)(|| 5 || 1 , 0O +7) 



hold on [0,T]. 

Proof: We use the transformation of variables 

Vh-yV = V(0,t,X,v) 



together with the representation for g from Lemma 2.2 and the estimate in 
Lemma |4.2j to obtain 



dv 



a(t,x)\ < / g(X(0,t,x,v),V(0,t,x,v)) 



H(V 2 (0,t,x,v))dv- / H(a 2 (t)v 2 )dv 



< 



g(X(0,t,x,v),V) Ideta^V^t^.u)! 1 dV 



H{V 2 )\detd v V(0,t,x,v)\ dV- a - 3 (t) H{V 2 )dV 



<C a - 3 (t)||5||oo+|^ 2 ) 



1 



\detd v V{0,t,x,v)\ a 3 (t) 



dV. 



Again by Lemma 4.2 we obtain the estimate 



1 



1 



\detd v V(0,t,x,v)\ a 3 (t) 



< 



1 



a 3 (t)\detd v V(0,t,x,v)\ 
•|detS„V(0,t,x,i;)-det(a(t)id)| 
< Ca~ 6 (t)Ca 2 (t)\d v V(0,t,x 1 v)-a(t)id\ 



<Ca- 3 (t) 1 , 

and inserting this into the above estimate for a we have 

IKi)||oo<Ca- 3 (i)(||g°||oc+7)- 

Since the £- integration extends only over the cube Q, the same estimate holds 
for ||cr(t)[[i. To obtain the estimate for d x a we differentiate the formula for g 
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in Lemma 2.2, integrate with respect to v, and use the same transformation of 
variables as above together with the estimates from Lemma 4.1 and 4.2: 



|0*<r(t,s)| <Ca(t) J (d {x , v) g)((X,V)(0,t,x,v)) 
+Cja(t) / \d v H(v 2 )\ v=V (o,t,x,v)\ dv 



dv 



<Ca- 2 (t)(\\g\\ hoo+1 ), 
and the proof is complete. 



□ 



5 The bootstrap argument 

The purpose of this section is to prove Theorem [0|. To do so we first note that 
under the assumption (FS7) we obtain an asymptotically stronger decay-in-time 
for the perturbed field than the one defining (FS7): 



Lemma 5.1 Let 7 be as in Lemma and let g be a .solution with g&T>, which 
satisfies (FS7) on some time interval [0,T]. Then the estimates 



and 



hold on [0,T] 



\\d 2 x W(t)\\oo<Ca- 13 / 16 (t) 



Proof: By Lemma 2.3 we have 

ll^w(i)iioo<c7o 2 (i)([k(t)iii /3 i|(7(i)||^ 3 +[k(t)ii 

and 



note that the right hand side of the Poisson equation ( l.lCj ) contains the factor 
a 2 (t). If we insert the estimates from Lemma 4.3 the assertion follows. 



□ 



Proof of Theorem |1 . Ij : Let e\ > and choose 

("1 



7< 



2d 



Next choose T > such that 



in such a way, that Lemma [4.2 holds; C\ denotes the constant from Lemma [4.3 
C 2 a 



(t)< ja- s (t), C 2 a- 13 / 16 (t)<la- 5 (t), t>T, 
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where C2 denotes the constant from Lemma 



5.1 



Finally, choose 



eo< 



ei 

2d 



in such a way that 



C 3 (T) £o + C 4 (T) e J 3/16 <^ a - 5 (T) 



where C3, C4 are the functions constructed in Lemma 3.1 and Lemma 3.3 re- 



spectively. Let g be a solution with initial condition g£T> such that ||<?||oo < £o- 
The functions C3 and C4 are nondecreasing and a is increasing. Therefore, by 
construction and Lemma [Tl] and |3.3| the free-streaming condition (FS7) with 
parameter 7 holds on some interval \0,T*\ with T* >T; we choose T* e]T,oc] 
maximal with this property. By Lemma ^l| the estimates 



and 



\\d x W{t)\\ 00 <C 2 a- 1 (t) 



\\d 2 x W(t)\\°o<C 2 a- 13 / 16 (t) 



hold on [0,T*[, and the choice of T implies that the free-streaming condition 
holds on [T,T*[ with parameter 7/2. Since T* is chosen maximal such that 
(FS7) holds on [0,T*[ it follows that T* = oo. In particular, 

lk(<)||oo<C 1 a" 3 (t)(||5°lloo+7)<eia" 3 W 



by Lemma 4.3 and by the choice of 7 and eg. It remains to see that this estimate 



implies the estimate on the density contrast S(t,x). Let (f,p,U) be the solution 
of the Vlasov-Poisson system J1.1Q, (|1.2[), ( |l~3| ) corresponding to (g,a,W), in 
particular, 

pit , a:) = po (t) + a(t , at) = p Q (t) ( 1 + a 3 (t)a(t , as)) . 

Thus 

6(t,x)=a 3 (t)a(t,x) 

where (g,a,W) is the perturbation written in the original variables, i. e., be- 



fore applying the transformation (1.8). If we denote the perturbation in the 



transformed variables by (g,a,W) then by fll.SD , 

a(t,x) = a(t,x) 

so that the estimate for <r, which we have established, implies the desired esti- 
mate for || <5(f) || oo; recall that we have been working in the transformed variables 
but have dropped the tildas for the sake of convenience. The proof of Theo- 
rem 1.1 is now complete. □ 



In addition to the assertion in Theorem 1.1 the following holds 
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Corollary 5.2 For every ei > there exists eo > such that every solution g of 



(|L9|), (1.10), (1.11) with g^V and ||g||i,oo<£o satisfies the estimate 

\\dMt)\\oo<a- 2 (t)ex, t>0. 



Proof: By Lemma |4 . 3| we obtain at the end of the proof of Theorem |l.l| the 
estimate 

||^a(i)||oo<C 1 a- 2 W(||. 9 || 1:OO + 7)<e 1 a- 2 (t), 

and the assertion follows. □ 

Final remarks: In the usual textbook explanations of the formation of galaxies 
it is argued that small perturbations of the homogeneous st ate lead to growth 
in time of the density contrast S(t,x) introduced in Theorem LI, cf. jn, 

Hill- 



However, we have proven that ||5(t)||oo can be made as small as we like, uni- 
formly in time, by making the initial perturbation small. This seeming con- 
tradiction can be resolved as follows: As already pointed out, our assumption 
a(0) = 1 means that we perturb the homogeneous state after is has already ex- 
panded out of the initial big bang singularity for some time. If we would pose 
initial data for the perturbation at times closer to the big bang then we would 
have to make these perturbations smaller in order to retain the same bound 
on ||<5(t)||oo. One reason for this is that if to < is defined by a(to) = then 
a(t) ~ (t — to) 2 / 3 for t >to close to to, regardless whether the energy E a of a is 
positive, zero, or negative. For E a = one has a(t) = C(t — to) 2 ^ 3 f° r a U t>to in 
which case the expansion would not be fast enough for our proof to go through. 

On the other hand, precisely this case E a = is considered in the "instabil- 
ity" arguments in the textbooks which for the generic case E a y^0 is justified 
under the restriction to times close to to, to which our analysis does not ap- 
ply. Mathematically, this initial growth of perturbations should thus not be 
referred to as an instability since it is not an assertion on the behaviour of the 
perturbations for all future time. 

After small perturbations have grown in the initial phase of the expansion the 
resulting solution, which on a larger scale can again be viewed as homogeneous, 
is now stable, provided we are in the case E a > 0, and no further conglomeration 
beyond the formation of galaxies and galaxy clusters occurs. 
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